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Abstract 

We study supersymmetric scenarios where the dark matter is the 
gaugino of an unbroken hidden U{1) which interacts with the visible 
world only via a small kinetic mixing with the hypercharge. Strong 
constraints on the parameter space can be derived from avoiding over- 
closure of the Universe and from requiring successful Big Bang Nucle- 
osynthesis and structure formation. We find that for typical values of 
the mixing parameter, scenarios with neutralino NLSP are excluded, 
while scenarios with slepton NLSP are allowed when the mixing pa- 
rameter lies in the range x ^ O(10~^'^ — 10^^"). We also show that if 
the gravitino is the LSP and the hidden U{1) gaugino the NLSP, the 
bounds on the reheating temperature from long lived charged MSSM 
relics can be considerably relaxed and we comment on the signatures 
of these scenarios at future colliders. Finally, we discuss the case of an 
anomalously small mixing, x 10^^^, where the neutralino becomes 
a decaying dark matter candidate, and derive constraints from gamma 
ray experiments. 



1 Introduction 



The existence of dark matter in the Universe is perhaps the most sohd indica- 
tion for physics beyond the Standard Model of particle physics [1]. Among 
the many extensions of the Standard Model that have been proposed in 
recent years, supersymmetry (SUSY) arguably remains as the most popu- 
lar. Among other merits, supersymmetric scenarios provide a very promis- 
ing candidate for the dark matter: the lightest supersymmetric particle 
(LSP) [21 [3]. With the particle content of the Minimal Supersymmetric 
Standard Model (MSSM), the LSP can be, over a large range of param- 
eters, either the lightest neutralino, the lightest sneutrino or the lightest 
stau. Among these, only the lightest neutralino is still allowed by present 
experiments as a viable dark matter candidate, provided i2-parity is almost 
exactly conserved. Furthermore, if supersymmetry is promoted to a local 
symmetry, the particle content of the MSSM also includes the gravitational 
supermultiplet, of which the spin 3/2 component, the gravitino, is also a 
viable dark matter candidate [1] , even if i2-parity is slightly violated [5l [6] . 

On the other hand, many extensions of the MSSM contemplate the possi- 
bility of a hidden sector, consisting of superfields which are singlets under the 
Standard Model gauge group. Hidden sector superfields usually couple very 
weakly to our observable sector, thus constituting a very natural arena for 
finding dark matter candidates. Generically, hidden sector particles couple 
to our observable sector only via non-renormalizable operators, presumably 
suppressed by a large mass scale, with a structure that is strongly model 
dependent. In consequence, deriving implications of the hidden sector dark 
matter on the thermal history of the observable Universe and for future col- 
lider experiments is hindered by our complete ignorance of the strength and 
the structure of the hidden sector interactions with our observable sector. 

There are however three instances where the hidden sector particles can 
couple to the MSSM particles via renormalizable operators, with a structure 
which is well defined by the Lorentz and gauge symmetries. Firstly, a hidden 
sector chiral superfield, S, could couple to the lepton and up-type Higgs 
doublet superfields via the Yukawa coupling SHuL in the superpotential, or 
to the two MSSM Higgs doublets, via SHuHd, provided these terms are also 
invariant under the hidden sector gauge group (as well as possible discrete 
and global symmetries of the theory). Secondly, if one of the MSSM chiral 
superfields is charged under a hidden sector gauge group, it will interact 
with the corresponding hidden sector gauge superfield, and in turn with 
other hidden sector chiral superfields via the i^-term. Finally, a hidden 
sector abelian vector superfield, X, may couple to the hypercharge vector 
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superfield through a kinetic mixing term, which is always allowed by the 
gauge symmetries [71 [8] . 

In this paper we will concentrate on the last situation and we will study 
scenarios where the hidden U{1) gauge group is unbroken at low energies. 
If this is the case, the corresponding hidden U{1) gaugino will have a mass 
comparable to the typical soft SUSY breaking masses of the observable sec- 
tor particles and, in some instances, smaller. Therefore, the hidden U{1) 
gaugino will not only be the lightest super symmetric particle of the hidden 
sector, but also the lightest among all the supersymmetric particles. 

Our motivation to consider an unbroken hidden sector U{1) group is 
twofold. First, in string theory compactifications hidden sector U{1) groups 
are ubiquitous and some of them could remain unbroken at low energies, 
in complete analogy to the familiar electromagnetic U{1) of our observable 
sector. Secondly, the case of the hidden sector unbroken U{1) is particu- 
larly intriguing, since this situation is practically unconstrained by present 
experiments. Indeed, it was shown long ago by Holdom that, in a non- 
supersymmetric world, the hidden U{1) gauge boson (the "paraphoton" [9]) 
completely decouples from the observable sector [10] . This result can be also 
generalized to a supersymmetric theory. Let us consider the SUSY invariant 
part of the Lagrangian, 

C = j (fe (wgWBa + W^Wxa + 2xWgWxa) + h.C. + 

+ j d^e<fe [<^^e^3YB^ ^ h^e^s^^h^ , (1) 

where the field strength superfield is defined as Wy = —jDDD°'V , V = 
B,X being the hypercharge or the hidden U{1) vector superfield, while 
$ and h denote, respectively, any Standard Model or hidden sector chiral 
superfield. Finally, x is the kinetic mixing parameter, which is induced 
through quantum effects by chiral superfields charged under both gauge 
groups. Without additional symmetries, values around x ~ 10"'^ — 10~^ are 
naturally obtained. However, e.g. in compactifications of heterotic |1H [T2] 
and type II [TSj [Ml [HI [TBI [17] strings, much smaller mixings are possible. A 
lower bound around x ^ 10~^^ was argued to hold in cases of gauge mediated 
supersymmetry breaking in heterotic string models [11], whereas in type II 
models with warped extra dimensions the kinetic mixing parameter can be 
exponentially small |17j . 

The gauge kinetic terms in Eq. ^ can be made canonical by introducing 
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shifted vector superfields, 



X = X + xB, 



B = ^Jl - X^B , 



(2) 
(3) 



leading to 




J dH {WgWBa + W^Wxa) + h.C. + 



(4) 



where g'y = Qy / \J^ — and g'-^ = XQx / ^/l — x^- Therefore, the canoni- 
cal normalization of the kinetic terms produces an unobservable shift of the 
hypercharge gauge coupling and the generation of a "minihyper charge" for 
the hidden sector chiral superfields [7]. Different astrophysical observations 
and laboratory experiments constrain the possible values of the minihy- 
percharge and the masses of the hidden sector particles. For instance, one 
obtains x ^ 10~^^ for masses below lOkeV from plasmon decay in red giants 
[181 119j (see also Ref. [2Q]). Nevertheless these bounds can be automatically 
avoided if the masses are large. In this case, as long as supersymmetry 
remains unbroken, the hypercharge vector superfield completely decouples 
from the observable sector and is not subject to any experimental constraint. 

The breaking of supersymmetry changes dramatically the previous pic- 
ture. Although the hidden U{1) gauge boson remains decoupled from the 
observable sector, we will show that in the presence of SUSY breaking effects 
a mixing between the hidden U (1) gaugino and the MSSM neutralinos is in- 
duced. Then, the unbroken hidden U{1) might produce observable effects 
in the cosmological evolution of the Universe, at collider experiments or in 
cosmic ray fluxes. 

Models with a hidden U{1) extension of the SM or MSSM have been 
extensively studied in the literature (see e.g. [21], for a recent review see 
Ref. [22]). Some of these works take into account kinetic- and mass-mixing 
[231 El 123 [261 [271 [281 [29], but it is typically assumed that the gauge symme- 
try of the additional U{1) is broken by a Higgs or Stueckelberg mechanism^ 
In this case, it is possible to derive bounds from high precision LEPI data 

'^Exceptions are e.g. Ref. [301 131) . where constraints on, and consequences of, higher 
dimensional operators that couple hidden and observable sector are studied. In Ref. [32] 
the authors study BBN and CMB constraints on the particle content of a completely 
decoupled hidden sector which may contain unbroken C/(l)s. Furthermore, see Ref. [33] 
for a short discussion about gauge coupling unification in the presence of kinetic mixing. 
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if the hidden U{1) gauge boson mass is large (namely, x ^ 0-05 for masses 
around 200 GeV [23] ) , and from different astronomical observations and lab- 
oratory experiments if the masses are small (namely, for masses around 
100 eV, the observed lifetime of the sun translates into a bound x ^ 10"^^ 
|34j . see Ref. |35] for bounds from considerations of the hidden CMB). 

The purpose of this paper is to consider in some detail cosmological 
constraints and phenomenological properties of a kinetically mixed hidden 
U{1) extension of the MSSM with unbroken gauge symmetry. We will as- 
sume that the hidden U (1) gaugino is the LSP in most of the paper, and 
concentrate on its prospect of being dark matter. In section 2, we will de- 
scribe the model in the component formalism and discuss typical values for 
mass- and mixing-parameters in scenarios with gauge and gravity media- 
tion of supersymmetry breaking. In section 3, we derive a bound on the 
mixing parameter from thermal overproduction of the hidden U{1) gaugino. 
Section 4 considers bounds from primordial nucleosynthesis for the cases of 
stau and neutralino NLSPs. There, we also analyze the effects of a possibly 
light gravitino and show how bounds on the reheating temperature can be 
relaxed. In section 5 we briefly discuss the collider phenomenology of this 
scenario. In section 6 we study the case with an anomalously small kinetic 
mixing, where the neutralino NLSP could become a viable, though unstable, 
dark matter candidate. Finally, we present our conclusions and an outlook 
in section 7. 



2 Model 

In this section, we will briefly describe the model under consideration in the 
more familiar component formalism. In the Wess-Zumino gauge, a vector 
superfield can be expanded in component fields as ^ = —Oa^OV^ + iOOOX — 
meX + ^eeeOD. Then, the pure gauge part of Eq. ([T]), including supersym- 
metry breaking effects, reads: 




where /C and Xi denote, respectively, the kinetic and mass mixing matrices: 

^) M=(f- f) . (6) 

\X ly V*^^ MbJ ^ ' 
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It is convenient to work in the basis where the kinetic terms are canonical. 
To this end, we will redefine the vector superfields according to 




(7) 



As explained in the introduction, the canonically normalized fields X^y and 
Dx completely decouple from the observable sector. On the other hand, 
a mixing between the hidden gaugino and the bino remains, unless 5M ~ 
xMx- More precisely, in the basis where the kinetic terms are canonical, 
the extended neutralino (5 x 5) mass matrix reads, to lowest order in x. 



(Mx 


SM 








\ 


6M 


Mb 





-Mzcpsw 


MzSf^sw 








Mw 


Mzcpcw 


-Mzspcw 





-Mzc/ssw 


Mzcpcw 







V 


Mzspsw 


-Mzspcw 


-/^ 


/ 



(8) 



where, 5M ~ 5M — x^x, Mx — Mx and Mb — Mb- Here, /x denotes 
the MSSM ^-term, Mz the mass of the gauge boson, sw the sine of the 
Weinberg angle and is related to the ratio of the two Higgs VEVs. 

Lastly, as was mentioned in the introduction, hidden sector matter par- 
ticles charged under the hidden U{1) will acquire a minihy per charge, as a 
consequence of the canonical normalization of the kinetic terms. For simplic- 
ity we will assume throughout the paper that hidden sector matter particles 
are too heavy to be produced during reheating, and thus do not infiuence 
the evolution of the Universe. 

The size of the gaugino mass terms Ai in the Lagrangian, Eq.([5]), is very 
model dependent. For example, in models with gravity mediation, gaugino 
masses would arise from the term 
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{zIWbWb + 2xZIWbWx + ZlWxWx). + h.c. , (9) 



Mp 



where coefficients of order one have been dropped. In this equation Zi are 
three spurion superfields, which will eventually acquire a vacuum expecta- 
tion value, (Zi) = Mp + O'^Fi, by some unspecified mechanisrrH, yielding 



^We assume that breaking of supersymmetry takes place in a sector that is distinct 
from the hidden 17(1) under consideration, hence the D-terms are zero and kinetic mixing 
has no effect on the MSSM mass spectrum (like e.g. in Ref. [11]). 
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Mb = Fi/Mp, 6M = xF2/Mp, Mx = Fs/Mp, where Mp ~ 2.4 x 10^^ GeV 
denotes the reduced Planck mass. 

The shnplest case consists on assuming just one spurion superfield which 
couples universally to all vector superfields, namely Zi = Z for all i. If this 
is the case, at the scale Mp, the different soft terms satisfy the relations 
Mb = Mx = ^-'^d 6M = \/3x"^3/2) '"^3/2 being the gravitino mass. 

Therefore, after the canonical normalization of the gaugino kinetic terms, 
the mass mixing term 5M in Eq. ^ vanishes. However, the universality of 
the gaugino mass terms is broken by quantum effects, thus inducing a non- 
vanishing 6M at low energies. Using the renormalization group equations 
for the gaugino mass matrix in appendix [Bl we find that the mass mixing 
at the electroweak scale reads 



where gx is the gauge coupling of the hidden U{1) gauge group, Mhid is the 
mass scale of hidden sector matter, and Bxx = J2i Qi sum, for all 

the particles in the theory, of all the hidden U (1) charges squared. On the 
other hand, the renormalization group evolution will also make the hidden 
gaugino mass smaller at low energies. If there is an enough number of 
(heavy) matter states in the hidden sector, the hidden gaugino will become 
the lightest super symmetric particle. 

In more general scenarios with gravity mediation, where several spurion 
fields contribute to the breaking of supersymmetry, a tree level mixing be- 
tween the bino and the hidden gaugino will usually remain after canonical 
normalization of the kinetic terms. On the other hand, the nature of the 
lightest supersymmetric particle, whether it is the gravitino, the bino or the 
hidden gaugino, depends on the details of the model. 

The most natural scenarios with a light hidden U{1) gaugino are proba- 
bly those with gauge mediated supersymmetry breaking. Assuming that the 
messenger fields are not charged under the hidden U{1) gauge group, the soft 
mass of the hidden gaugino will entirely come from gravity effects. If this 
is the case, the gravitino and the hidden ^7(1) gaugino will acquire similar 
masses, which are naturally much smaller than the masses of the supersym- 
metric particles of the MSSM. Again, whether the LSP is the gravitino or 
the hidden gaugino depends on the particle content of the hidden sector, 
which will drive the hidden gaugino mass to smaller values at low energies 
through radiative effects. 

An important quantity in our analysis will be the mixing angle O between 




(10) 
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the bino and the hidden U{1) gaugino mass and interaction eigenstatei 
The approximate diagonahzation of Eq. ^ yields 



where we have emphasized that all the quantities should be evaluated at the 
electroweak scale. In the scenarios of supersymmetry breaking proposed in 
this section with hidden U{1) gaugino as LSP, the mixing angle Q can be 
written as 

e-c.xEw-g|wf^. (12) 

where C is a constant that can roughly vary between 10"^ < C < 1, de- 
pending on whether (5M^^ is already present at tree level or is generated 
radiatively, and on the particle content of the hidden sector. For definiteness 
we will choose C = 1 in the subsequent analysis. 

Following the above discussion, we will first study the phenomenology of 
a hidden U (1) gaugino LSP, and second the case where both, the hidden U (1) 
gaugino and the gravitino, constitute the lightest supersymmetric particles. 



3 Overproduction constraints 

The least model dependent bound on the hidden U{1) gaugino scenario 
comes from considerations of the thermal history of the universe. After 
reheating of the visible sector, hidden U (1) gauginos can be produced in the 
hot primeval plasma through the mixing with the MSSM neutralinos. 

This mechanism is similar to the thermal production of gravitinos or ax- 
ions (see e.g. Ref. [36l EH [38] ) . However, in contrast to these cases, the final 
abundance of the hidden U{1) gaugino does not depend on the reheating 
temperature, Tr. This follows from the fact that mass- and kinetic-mixings 
appear at the renormalizable level of the Lagrangian, whereas gravitino and 
axion couplings are mass-suppressed. Actually, the production is most effi- 
cient at low temperatures, and this makes an exact calculation of the hid- 
den U{1) gaugino abundance extremely complicated, since at temperatures 
around T 100 GeV the QCD coupling constant gs is of the order of one, 
which precludes a sharp separation between hard, T, and soft, gsT, momenta 
as required for the proper treatment of hard thermal loops and screening 

^Note that we always assume the mass hierarchy Afs < Mw < fi in this paper. 
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Process 


\M\ye-'{g'gsr\nj-Yl^ 


1) qg qXx 


2— 


2) qg qXx 


-2(! + f) 


3) qg qXx 




4) qg qXx 


4 


5) qq gXx 




6) qq gXx 


8 


7) qq gXx 





Table 1: Squared matrix elements for inelastic two-body scatterings with a 
hidden U{1) gaugino in the final state. Here, T^^ denotes the generators of 
the fundamental representation of the SU{3) gauge group, 11,/? the hyper- 
charges of the (s)quarks, and g' and gs the gauge couplings of the electroweak 
and strong interactions, respectively. We average over initial and sum over 
final spins. 



effects ^9]. How to deal with this drawback is an open problem even in the 
milder case of gravitino production (see, however, Ref. |38|). 

For the purposes of this paper it is enough to derive an order of mag- 
nitude estimate on the hidden U{1) gaugino abundance. To this end, we 
have calculated the collision integral in a relativistic QCD plasma, taking 
into account 2 — > 2 scattering in the QCD and hypercharge sector with a 
bino in the final state. Since binos mix with hidden U{1) gauginos with 
the mixing angle G, the production rate of hidden f/(l) gauginos can be 
straightforwardly calculated by multiplying the result for the binos by 0^. 
The relevant scattering processes and the corresponding squared matrix el- 
ements are listed in Tab. [H 

Summing over all spins and particles, the collision term of hidden U{1) 
gauginos in the relativistic MSSM plasma read^ 

~ 3 X io"^e(T)2r^ , (13) 

where one has in general to take into account the temperature dependence 
of the mixing angle 0, which stems from the thermal mass of the bino (see 
Eq. (I12p ). Here, we simply take Mb{T) ~ Mb, since most of the production 
is expected to happen when the particles become non-relativistic. In this 

*We use Boltzmann distribution functions for simplicity. 
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equation, a' and as denote the hypercharge and QCD couplings, respec- 
tively, and jE — 0.577. Furthermore, k* denotes the cut-off of intermediate 
three-momenta in t-channel processes. Only processes with quarks in the 
intermediate state contribute to these divergences, and we identify k* with 
their thermal mass (at T ~ 100 GeV the thermal quark masses lie between 
0.63r and 0.84T [38]). Finally, solving the Boltzmann equation yields the 
thermal abundance of hidden ^7(1) gauginos. 



Taking into account thermal masses could actually lead to a further order- 
one enhancement due to particle decay in the QCD plasma [38]. Contribu- 
tions from electroweak interactions and Yukawa couplings are also expected 
to give sizable corrections, again of order one. 

Requiring that Q.xh^ ^ 0.1 gives the overproduction bound 



where we have used Tq ~ Mb as a cutoff in Eq. (|14p . which corresponds to 
squark- and gluino masses around 2>A4b- For other values the bound scales 
roughly like \/To. 

In Fig. m (for Mb = 180 GeV) and Fig. [2] (for Mb = 150 GeV) we show 
the region in the (x,Mx) parameter space where the hidden U{1) gaugino 
is overproduced. Remarkably, a large part of the parameter space suggested 
by string theory is excluded by the constraint Eq. (fT5|l H 

4 Bounds from primordial nucleosynthesis and 
structure formation 

In general, dark matter can be produced thermally and/or in the late decay 
of some other relic particle ^. The abundance of dark matter today is then 
given by 0dm = ^^dm + {^'Dm/M^)^1^^- , where refers to the thermal 
abundance of the particle ^. Typical candidates for ^, like the stau or a 
neutralino, naturally freeze out with an abundance of the right order of 
magnitude to account for the observed dark matter abundance, and a given 

sizable annihilation rate between hidden C/(l) gauginos would allow to circumvent 
this bound. If the mixing angle Q is small, this would require additional light matter 
states in the hidden sector, which do not exist in our scenario. 




(14) 




(15) 
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ii) 


Hi) 


iv) 


NNLSP 






slepton 


slepton 


NLSP 


slepton 


neutralino 


gravitino 


Ax 


LSP 


Ax 


Ax 


Ax 


gravitino 



Table 2: The different scenarios that we will investigate in this work. 

dark matter candidate can inherit this property as long as its direct thermal 
production is negligible and provided the masses of the two particles are 
comparable. 

These so-called super WIMP scenarios (super- weakly interacting massive 
particles, see Ref. [M]), where the dark matter super WIMPs are produced in 
the late decay of other particles, are subject to strong constraints due to their 
potential impact on primordial nucleosynthesis and structure formation. If 
the production takes place during or after Big Bang Nucleosynthesis (BBN) , 
the hadronic and electromagnetic energy released in the decay can destroy 
the successful predictions of the standard scenario (see e.g. Ref. [501133]). In 
cases where the decaying particle is electrically charged, it could additionally 
form bound states with "^He, triggering the catalytic production of ^Li and 
^Be |4H 144] , Furthermore, the superWIMP itself is potentially produced 
with a momentum that is large enough to wash out density fluctuations 
on small scales with impact on structure formation, similar to scenarios 
with warm dark matter (for a discussion in the context of superWIMPs see 
Ref. [45]). 

Below, we will apply these bounds to the case of a hidden U{1) gaugino 
superWIMP, concentrating on the four different scenarios shown in Tab. [21 
Following the superWIMP paradigm, the relic abundance of the neutralino 
or the slepton is fixed by the requirement that the LSP constitutes all dark 
matter today. 

Scenario i) After freeze out, the slepton NLSP decays predominantly 
through the channel / — > Ax^ into a lepton and the hidden U{1) gaugino 
LSP. The corresponding decay width is given by 

r- ~ ^e^y-^Af- ( i - ^\ riei 

where and Mj denote, respectively, the hypercharge and the mass of the 
slepton NLSP. 
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The decay induces electromagnetic cascades in the plasma,_ the branching 
fraction into electromagnetic energy being close to one, B^^y^ ~ 1. The 
corresponding energy release is given by e^j^ w {M? — M^) / {2Mj) (see 
Ref. [l6] for an exhaustive discussion). 

Contributions to the hadronic energy release stem, if kinematically al- 
lowed, from three-body decays into Z^, and Higgs bosons. The corre- 
sponding branching ratios are all of the order of, or smaller than, a' / (IGvrs^) ~ 
10~^, whereas the hadronic branching fractions of the subsequent decay pro- 
cesses are all of the order of one. If these decay channels are kinematically 
closed, the dominant contribution to hadronic energy release comes from 
the strongly suppressed four-body decay / Xxlqq, with a branching ra- 
tio of the order of 10~^. Hence, the overall hadronic branching fraction is 
-^had ~ 10"^ foi' large masses of the hidden U{1) gaugino {Mx ^ M^- — Mz) 
and B!^^^ ~ 10~^ for small masses {Mx ^ Mj — Mz)- Variations in the 
corresponding energy release e\^.^^ are subdominant because i^^ad varies by 
several orders of magnitude, and we simply take e^^^^ = (1/3) (M^ — Mx) 



Applying the BBN bounds from Ref. |l0] to this scenario, a consider- 
able part of the hidden U{1) gaugino parameter space can be excluded (see 
Fig. (Tjlfl We also show the region that is excluded by the catalytic enhance- 
ment of ^Li and ^Be production. These effects essentially require that the 
lifetime of the slepton is smaller than around 2 x 10^ s |41j . 

The impact of super WIMP scenarios on structure formation is encoded 
in the free streaming length Afsi defined as the comoving distance traveled 
by a particle since it was produced: 



On scales below the free streaming length of dark matter particles, den- 
sity fluctuations are washed out before structure formation begins (see e.g. 
Ref. [IZ])- The non-observation of this effect in the latest Lyman-a Forest 
data from the Sloan Digital Sky Survey [48] excludes a free streaming length 
larger than roughly Aps ^ O.5Mpc0 The region in the parameter space of 
the hidden U{1) gaugino that is excluded by this requirement is shown in 
Fig. [Hand lies completely in the region already ruled out by BBN. 

''In Ref. |40) the bounds are shown for the cases Bhad = 1, 10"^ and zero. We just 
interpolate between Shad ~ 1 and Bhad ~ 10"'^ hnearly in logj^Q(_Bhad) and use B^ad = 
as a cutoff. 

'^See Ref. [49], where the power spectrum of non-thermally produced dark matter is 
confronted with lower bounds on the mass of sterile neutrino dark matter |50l 151) . 



m 




(17) 
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Scenario ii) Two-body decays of a neutralino NLSP lead to final states 
where the hidden U{1) gaugino LSP is accompanied either by a or Higgs 
boson or a photon. However, if the decay into or Higgs bosons is kinemat- 
ically forbidden, neutralinos predominantly decay via — > ^xff, since the 
radiative two-body decay into photons is one-loop suppressed and negligible 
in most cases. We have summarized all relevant decay widths in appendix lAl 
particularized to the case of a bino-like neutralino. The hadronic branch- 
ing fraction for the decay of a bino-like neutralino is of the order of one 
for hidden U (1) gaugino masses that allow the decay into bosons, and 
can range between 10~^ and order one otherwise (depending on the squark 
masses). For definiteness we choose B^^^ ~ 0(1) for all masses Mx, since 
smaller values at large Mx will not change our conclusions. For the corre- 
sponding energy release we simply take ehad = 1/3(Mb — Mx)- We apply 
the BEN bounds of Ref. [30] and show the corresponding excluded region 
in Fig. O Together with the overproduction bounds, this scenario is ex- 
cluded in the whole {x,Mx) parameter space suggested by string theory, 

Scenario iii) We will now discuss the effects of an additional light grav- 
itino, concentrating on the spectra shown in Fig. [3l where we also show 
typical values for the relevant decay widths. 

If the hidden U{1) gaugino is the LSP (see Fig. [3ti), it can be produced 
in three different ways: thermally, via slepton decay and via gravitino decay. 
For mixings and masses that are allowed in Fig. [H the decay of sleptons into 
gravitinos is strongly suppressed, and the decay of gravitinos into hidden 
U{1) gauginos is only accompanied by hidden U{1) gauge bosons and invis- 
ible to MSSM particles. Hence, bounds from BBN and overproduction are 
essentially the same than in case i). 

However, this does not hold for the free streaming bounds, since, in con- 
trast to hidden U{1) gauginos which are produced thermally or via slepton 
decay, hidden U{1) gauginos that stem from the decay of gravitinos can have 
a quite large free streaming length, of the order of several Mpc, as shown 
in Fig. m This leads to dark matter with a cold and a warm component, 
so-called ACWDM models. 

Constraints on ACWDM models where recently studied in Refs. \52\ [531 
I54j . A warm dark matter component would induce a step in the power 
spectrum of density fluctuations, with a size which depends on the fraction 
of the dark matter that is warm and a position which is roughly given by 
the corresponding free streaming length as /cps ~ 27r/AFs. Although models 
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with mixed cold/warm dark matter arise quite naturally in many situations, 
a full general analysis with all the latest data is still lacking|f| Bearing in 
mind this limitation, we will simply show the corresponding bounds on the 
reheating temperature and mixing parameter for the two reference values 
/ = 0.2 and / = 0.02, where / denotes the fraction of dark matter that is 
warm, r^oM = /f^WDM + (1 — f)^CDM] a complete analysis of the mixed 
cold/warm dark matter scenario is beyond the scope of this paper. 

As discussed above, the warm component of hidden U{1) gaugino DM 
comes from the late decay of gravitinos. The thermal abundance of graviti- 
nos directly depends on the reheating temperature Tr like |36^ [371 EH] 

where mg denotes the gluino mass. Hence, the bounds on ACWDM models 
translate into a bound on the reheating temperature as shown in the right 
part of Fig. [5l 

Scenario iv) Scenarios with a gravitino LSP and a slepton (usually a stau) 
NLSP are known to be strongly constrained due to the catalytic production 
of ^Li and ^Be in the presence of long-lived charged particles during BBN (see 
above discussion). The resulting bound on the lifetime of staus translates 
into an upper bound on the gravitino mass of around ttiq < 1 GeV. At the 
same time, requiring that gravitinos with ttIq < 1 GeV are not overproduced 
forbids reheating temperatures as high as Tr ~ 10^ GeV (see left part of 
Fig.[5]and Eq. ()18p ). However, these high reheating temperatures are favored 
in scenarios with leptogenesis as the source of the baryon asymmetry of the 
universe [561 EZ] • 

One way to relax the tension between leptogenesis and gravitino DM 
is to assume that the production of entropy between stau freeze out and 
BBN dilutes the stau abundance sufficiently to evade the bounds |581 159] . 

*In Ref. [54] a consideration of the uncertainty of the power spectrum from WMAP 
1-yr data [55] bounds the fraction / of dark matter that is aUowed to have Afs ~ 6 Mpc to 
be / < 0.2 — 0.4. This seems to be consistent with a hydrodynamical analysis in Ref. [52j . 
where the bound on the fraction of dark matter that can be made out of thermal light 
gravitinos with a mass of tUq « 16 eV is stated as / < 0.12, since the corresponding free 
streaming length (using Eq. |T7J)) is around Afs ~ 40 Mpc. However, both analyses do 
not take into account the latest Lyman-a-forest data [48]. An analysis in Ref. [53] which 
incorporates these data, and where lower bounds on the mass of sterile neutrino dark 
matter from Ref. [50) are rescaled for the case of ACWDM models, suggests that even a 
component with a free streaming lengths around Aps ~ 4 Mpc is bounded strongly with 

/<o.i. 
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Another way is to impose a small violation of i2-parity [6] , which allows the 
stau to decay fast into standard model particles. In this case, the gravitino 
becomes unstable and decays with a lifetime which is typically longer than 
the age of the universe. The gravitino decay products could be detected as 
an anomalous contribution to the cosmic ray fluxes [001 ED E2] , opening the 
possibility to probe this scenario. 

Interestingly, a similar mechanism can work in the presence of a hidden 
U{1) gaugino for spectra of the form (slepton. Ax, gravitino), as depicted 
in Fig. [3b. For mixings and masses of the hidden U{1) gaugino that lie in 
the allowed region of Fig. [H the lifetime of the stau is compatible with all 
BBN bounds. In this case, the reheating temperature can be as large as 
Tr ~ 10^ GeV, as shown in Fig. [5l since the gravitino mass and the stau 
lifetime are decoupled. 

The warm component of gravitino dark matter, with a free streaming 
length of the order of several Mpc (see Fig.[4|), comes now from the late decay 
of hidden U{1) gauginos. Hence, bounds on ACWDM models constrain the 
sum over the thermal and non-thermal abundance of hidden U(l) gauginos 
in this case. Firstly, this shrinks the allowed region of the hidden U{1) 
gaugino parameter space as shown in Fig. [T](dashed lines). Secondly, in cases 
where the thermal production of hidden U{1) gauginos is negligible, the stau 
abundance is bounded like Fstau < (9 x W''^^ - 9 x W-^^){M^ / WO GeV)~'^ 
if we allow a fraction / = 0.02 — 0.2 of dark matter to be warm. Thus, the 
^Be and ®Li bounds on the abundance of long-lived staus, which are roughly 
given by Kjtau < 10~^^ and Kjtau < 10"^^, respectively (see Ref. [44]), can 
be weakened by two or three orders of magnitude. 

In the CMSSM with gravitino dark matter the typical relic abundance 
of stau NLSPs ranges between l^tau ~ 10~^^ and ygtau ~ 10~^^ (although 
considerably smaller values are possible due to resonant effects or coanni- 
hilation [63], or e.g. in models with non-universal Higgs masses [M])- In 
these kind of models, the phenomenologically allowed parameter space will 
be greatly enlarged by the above mechanism. However, a detailed analysis 
of the influence of a weakly interacting hidden sector on constraints on the 
reheating temperature would require reliable bounds for general ACWDM 
models. 

5 Collider phenomenology 

At the LHC, cascade decays of squarks and gluinos can produce of the 
order of 10^ NLSPs per year if the sparticle masses are close to the current 
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experimental limits [65] , thus offering a miique opportunity to test scenarios 
with a light hidden U{1) gaugino. It was shown in the previous sections that 
in this class of scenarios the neutralino as NLSP is excluded, whereas the stau 
as NLSP is allowed only if it has a lifetime that ranges between r ~ 10~^ s 
and T ^ 10^ s. Therefore, the stau will typically traverse the detector leaving 
a distinctive heavily ionizing charged track [66]. Moreover, if the stau is slow 
enough, it might get trapped in the detector, decaying at late times and thus 
allowing a measurement of the lifetime and, correspondingly, of the mixing 
parameter x- The number of staus trapped can be incremented by placing 
a 1-10 kton massive material around the LHC detectors, which will allow to 
collect 0(10^ — 10^) staus [67j. If the stopper material is simultaneously an 
active real-time detector, the decay products and their angular distribution 
could be studied in detail. This will allow to distinguish this scenario from 
the scenario of gravitino dark matter, which yields very similar signatures 
at colliders [68l . 



6 The case with an anomalously small x 

We have shown in sections [3] and H] that constraints from overproduction, 
BBN and structure formation exclude the scenario with a neutralino NLSP 
and hidden U{1) gaugino for natural values of the mixing parameter x ~ 
(see Fig. [2]). However, these strong constraints could be evaded 
in more elaborated models yielding a smaller x- Our interest in exploring the 
region with small x is that the neutralino lifetime could become larger than 
the age of the Universe, thus constituting a viable dark matter candidate 
itselfo Nevertheless, in contrast to the standard neutralino dark matter 
scenario, in this scenario the neutralino NLSP is unstable and decays with 
very long lifetimes into the hidden gaugino and Standard Model particles. 

The neutralino dark matter particle (that for definiteness we assume to 
be a pure bino) can decay into different channels: 

Xb Z Xx , 

Xb 1 , 
Xb h Xx , 

Xb ^ //Ax, (19) 

® Clearly, the scenario with charged slepton NLSP is excluded for very small values of 
X, since it would yield an abundance of anomalously heavy hydrogen in conflict with the 
experimental constraints. 
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with decay rates listed in appendix |Al The fragmentation of the and 
the Higgs boson produces a continuous spectrum of stable particles, such as 
photons, positrons or antiprotons which contribute to the primary fluxes of 
cosmic rays. The experimental constraints on these fluxes will then translate 
into constraints on the parameters of the model. 

On the other hand, it has been pointed out recently that the decay of a 
dark matter particle with a lifetime r ~ lO^^s into weak gauge bosons with 
a momentum ~ 50 GeV could simultaneously explain the EGRET anomaly 
in the extragalactic gamma-ray background and the HEAT excess in the 
positron fraction |6HI62j. Interestingly, the scenario we are considering here, 
where the neutralino dark matter particle decays into a gauge boson, 
constitutes a promising candidate of this kind. Moreover, in contrast to 
other candidates for decaying dark matter that have been proposed, such as 
the gravitino [6j or a hidden U{1) gauge boson [69], the unstable neutralino 
could be detected in direct dark matter searches. 

Although the simplest models with hidden gaugino LSP and neutralino 
NLSP predict a mixing parameter x ranging between 10"'^ and 10~^^, and 
thus a neutralino lifetime ranging between r ~ 10"^^ s and lO^S, it IS 
not difficult to construct plausible models where the neutralino lifetime can 
be orders of magnitude longer than the age of the Universe. To this end, let 
us consider a model with two hidden U{1) gauge groups, U{l)x and U{l)x', 
and vector superfields denoted by X and X' , respectively. We will assume 
that the hidden sector particles that generate the kinetic mixing are charged 
either under U{1)y and U{l)x', or under U{l)x and U{l)x', but not under 
U{1)y and U{l)x simultaneously. If this is the case, the kinetic mixings 
W§Wx'a and WxWx'a wiU be generated, but not W^Wxa- Therefore, 
the SUSY invariant part of the Lagrangian reads 



-^SUSY 



/ 



cf0 {wgWBa + W^Wxa + W^'Wx' a + 




(20) 



In addition, there exists a soft SUSY breaking Lagrangian: 




(21) 



where 




(22) 



17 



This structure for the soft mass matrix can be straightforwardly derived 
from Eq. (j20p introducing five arbitrary spurion superfields as in Eq. ([9]). 
Note that JMi and 5M2 are suppressed by xi ™d X2-, respectively. 
In order to canonically normalize the vector superfields, we define 

B = B , 

X = X-X1X2B, 

X' = X' + xiB + X2X, (23) 

where we have assumed X1tX2 ^ 1- Then, the SUSY invariant part of the 
Lagrangian reads: 

CSVSY = J {WgWsa + W^Wxa + W^,Wx' a) + h.C. , (24) 

while the SUSY breaking Lagrangian reads: 

1 f^A 
CsmY=--{Xx \x' Ab)A^ I Ax' I +h.c. , (25) 

with 

/ Mx 6M2 5Mi2\ 
M ~ 5M2 Mx' 5Mi . (26) 
\SM12 SMi Mb I 

In this equation, Mx — Mx, Mx' ~ Mx' and Mb — Mb, while 6M2 ^ 
5M2 — X2MX' and 5 Ml ~ 6M1 — xi^x- More importantly, after the 
canonical normalization of the kinetic terms, a mixing term between the 
bino and the hidden gaugino Ax has been generated, 6M12 ^ Xi X2 {Mx + 
Mx') — X2 SMi — xi 6M2, which is doubly suppressed by xi and X2- 

If Ax is the LSP and A^ is the NLSP, the decay rate of the bino NLSP is 
very strongly suppressed, thus yielding very long lifetimes for the neutralino 
dark matter. For instance, assuming that the kinetic mixing between the 
observable sector and the hidden sector is xi ~ 10~^^, as suggested by string 
theory, a mixing between the two hidden sector U{l)s of X2 ^ 10"'^ would 
be necessary in order to render a neutralino lifetime of the order of 10^^ s. 

The signatures in cosmic ray experiments of the scenario of neutralino 
dark matter which decays into hidden gauginos are fairly model dependent 
and will be shown elsewhere [70]. In this paper we will just derive a conser- 
vative bound on the mixing parameter x from requiring a gamma ray flux 
from neutralino decay in agreement with the EGRET observations. To this 
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end, we will just consider the contribution to the gamma ray flux from the 
two body decay — > 7Ax- 

The total gamma ray flux receives several contributions. First, we expect 
a gamma ray flux from the center of the Galaxy produced by neutralino 
annihilations. This contribution does not depend on whether the neutralino 
is absolutely stable or on the contrary decays at late times. Hence, it does 
not constrain our scenario and will not be further discussed. The interested 
reader can find in the vast existing literature the prospects of detecting such 
a signal from neutralino annihilation coming from the center of the Galaxy 
(see e.g. Ref. [1]). 

In addition, we expect a diffuse flux produced in the decay of neutralinos 
at cosmological distances and in the Milky Way halo. Following Ref. [6l I71j. 
we will average the halo gamma ray signal of the — > 7 Ax channel over 
the whole sky, excluding a band of ±10° around the Galactic disk in order 
to compare our results with the EGRET data [72]. To be conservative, 
we will require that the corresponding peak in the energy spectrum, after 
convolving with the energy resolution of 15% quoted by EGRET, remains 
below the 2a band of the EGRET spectrum obtained by Sreekumar et al. 
in [73], 

T - f"^-"-^^) - '»-'<-^ ^ (045^)"""°' ■ 

(27) 

The result is shown in Fig. [6l where we used a reference bino mass of Mb = 
150 GeV. For hidden U{1) gaugino masses around Mx ~ 75GeV, mixings 
down to X ~ 10~^^ can be excluded from the EGRET measurements of the 
gamma ray flux. 



7 Conclusions 

An unbroken hidden U{1) that interacts with the Standard Model only via 
kinetic mixing with hypercharge decouples completely from the visible world. 
However, in the supersymmetric version of this scenario this is no longer the 
case. We have shown that a mass mixing between the bino and the hidden 
U{1) gaugino is always generated via radiative effects, although this mixing 
can be generated already at tree level in some well motivated scenarios. 
Moreover, we have discussed different scenarios of supersymmetry breaking 
in which the hidden U{1) gaugino is the lightest supersymmetric particle. 

We have mostly concentrated on this possibility and we have derived cos- 
mological bounds on this scenario from precluding overproduction of hidden 
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U{1) gauginos and from the requirements of successful Big Bang nucleosyn- 
thesis and structure formation. The combination of these constraints ex- 
cludes a neutralino NLSP except for extremely small mixings (see Fig. [2]). 
On the other hand, when the NLSP is a stau, an allowed window for mixings 
around x ~ 10~^^ — 10"^^ remains (see Fig. [T]). In this window, the stau has 
a lifetime larger than 10~^ s and thus might be detected at future colliders 
as a heavily ionizing charged track. 

The reheating temperature in scenarios with gravitino dark matter and 
stau NLSPs is known to be strongly constrained, and we have shown that 
these constraints relax considerably in the presence of a hidden U{1) gaugino 
(see Fig. [5]). This might be a rather general effect of very weakly interacting 
hidden sectors and deserves further attention. 

Finally, we have discussed the case of an anomalously small mixing pa- 
rameter, X 10~^^. For these small mixings, the neutralino NLSP can 
become long lived enough to constitute the dark matter of the Universe. 
We have constructed a simple model with two hidden U{l)s where a tiny 
mixing can be naturally obtained. Even though the neutralino is very long 
lived, it eventually decays into the hidden gaugino and standard model par- 
ticles, which might be detected as an anomalous contribution to the cosmic 
ray fluxes. Using the EGRET measurement of the extragalactic gamma 
ray flux, we have derived a conservative bound on the mixing parameter 
X < 10~20-2i (ggg Fig. ED. 
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A Decay rates 

For convenience, we summarize the neutralino, slepton and gravitino decay 
rates that we used in this paper. 

The decay widths for bino-like neutralinos that decay into hidden U{1) 
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gauginos can be derived from the results in Ref. |74l I75j . They are 

V327r2 J 
t ( ^ 



'B 

M% (Ml Ml (Ml MIS' MxMb 

fi^Ml ^ Af| ^ M| ^ \Ml M| y ^ Mz Mz 

XlT^I I-ITJT 1-T# > (29) 



1 ,0^0,. MUL,,fo Ml Ml 

' Ml' Ml 



nXB^hXx) - -k^9"Q'MB^^X'l'{l 



Mi M^ Mx\ 

V{\B^l!\x) ^ L4xlO-V*e2AfB||(l-^J)' , (31) 

where we have used the function 

X{a^,b^,c^) = (a^ - {h + cf){a^ - {h - cf) . 

The approximations that we use for the neutralino mixing angles become ex- 
act in the hmit Mx ^ Mb Mw <C ^. Note that the two-body decay into 
a hidden U{1) gaugino and a photon, Eq. (p9]) . is one-loop suppressed, and 
we sum over ah (s)leptons in the loop, assuming that they have a common 
mass Msf. The final state in Eq. (jSip incorporates neutrinos and charged 
leptons, and we took into account sfermions and gauge bosons in the 
intermediate state, where the later are subdominant. 

The relevant decay widths with gravitinos in the initial or final state can 
be found in Ref. [461 [76] and are 

T(a^XxX) = -J-jM? I 1 - 11^. I + (33) 
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B Renormalization group equations 



The RGEs for multiple U{1) models with kinetic mixing where first pub- 
lished in Ref. [77]. One-loop RGEs for gaugino masses and their mixings 
can be found in Ref. [78]!^ 

The RGEs acquire their simplest form in the basis where gauge bosons 
and gauginos have a canonical kinetic term. The matrix of the gauge cou- 
plings Qij is defined according to the term C D gijj^Ali in the Lagrangian, 
where the indices i = 1,2 and j = 1,2 run over the charged currents of the 
two sectors and over the gauge boson states, respectively. In this represen- 
tation, the one-loop RGEs for the coupling constants are 

jf9 = jt^afB-g , (35) 

where we have used the charge matrix Bij = tr{QiQj), and t = ln(Q/(5o) 
with Q as RG scale. The trace in Bij runs over all chiral supermultiplets. 
The RGEs for the gaugino mass matrix A4 take a similar form 

Jt^^lt^ ^^^''Bg + g'^BgM) . (36) 

Note that only three of the four entries of g are physical because the basis 
of the gauge bosons is only fixed up to a rotation. It is convenient to state 
the RGEs for the case where the non-diagonal elements in the couplings gij 
are chosen to vanish, and where their third free component is absorbed in 
the kinetic mixing term like in Eq. ([5]). We obtain 

Jt9x,B = ^9\,bBxx,bb , (37) 

= -^gxgsBxB + {gxBxx + qIBbb) + O(x^) , (38) 

^Mx,B = -^gl,BBxx,BBMx,B , (39) 

^<5M = {glBxx + qIBbb) SM + 0{x^) . (40) 
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Note that our result for the RGEs of the gaugino mass matrix in Eq. p6|) differs from 
the one given in Ref. [78] by an additional symmetrization. 
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Figure 1: Summary of bounds on the hidden U{1) gaugino parameter space 
for the case of a slepton NLSP. We use = 150 GeV and Mb = 180 GeV. 
The upper dark blue region is excluded by thermal overproduction. Below 
this region, the hidden U{1) gaugino is dominantly produced via late decay- 
ing sleptons. The light blue region is excluded by energy injection during 
BBN |40] . whereas the light green region is excluded by catalysis of ^Li pro- 
duction [H]. We also show the region which would be excluded solely by 
free streaming arguments. The dotted lines show the slepton lifetime. In 
the presence of a gravitino with Mq = 100 GeV the slepton would domi- 
nantly decay into the hidden U{1) gaugino, except in the red lower region. 
The dashed lines show the region that is potentially excluded by bounds on 
AC WDM models in scenarios where the thermally produced hidden U{1) 
gaugino decays into a gravitino LSP with large free streaming length (see 
Fig. HI). The lines correspond to a dark matter fraction / = 0.02 and / = 0.2 
which is allowed to be warm. 
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Figure 2: Summary of bounds on the hidden U{1) gaugino parameter space 
for the case of a bino-hke neutrahno NLSP. We take Mb = 150 GeV. The 
dark blue region is excluded by thermal overproduction. Below this region, 
the hidden U{1) gaugino would be dominantly produced via late decaying 
neutralinos. This scenario is totally excluded by BBN [ID] (light blue re- 
gion). The bound actually strongly overlaps with the overproduction region 
(dashed line). We also show the region (in yellow) that would be solely ex- 
cluded by free streaming arguments. The dotted lines show the lifetime of 
the neutralino. (We used Mgf = 400 GeV and fi = 300 GeV for the branching 
ratios) . 
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Figure 3: Spectra for (a) gravitino NLSP and (b) gravitino LSP (c/. last 
two cases in Tab. [2|). We show the widths for the different decay processes 
for typical particle masses. The mixing parameter is assumed to lie in the 
allowed region of Fig. [TJ 




Figure 4: Contour plot of free streaming lengths in units of Mpc. The upper 
left (lower right) corner shows the free streaming length of gravitinos (hid- 
den U{1) gauginos) that stem from the late decay of hidden U{1) gauginos 
(gravitinos) . 
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Figure 5: Bounds on the reheating temperature as function of the gravitino 
mass, using Eq. (jlSp with nig = 800 GeV. The mass of the hidden U{1) 
gaugino is fixed to Mx = 120 GeV. If the gravitino is the LSP (left part), 
the reheating temperature is only bounded by overproduction arguments 
(dark red region), which are only slightly strengthened when, say, 20% of the 
gravitino abundance is due to non-thermal production (light red region). A 
gravitino NLSP (right part) would late decay into the hidden U{1) gaugino, 
yielding a warm dark matter component. If only a fraction of 20% or 2% of 
dark matter is allowed to be warm (with free streaming lengths as shown in 
Fig. H]), the corresponding blue regions are excluded. 
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Figure 6: Region of the hidden U{1) gaugino parameter space which is ex- 
cluded by confronting the EGRET measurements [73] with the extragalactic 
gamma ray flux produced in the decay Ab — > jXx (blue region). The bino 
mass is fixed to Mb = 150 GeV. We also show for reference the isocurves 
for the neutralino life-time. 
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